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A Semi-Quantitative  Model  for  the  Burning 
Rate  of  Solid  Materials 


James  G.  Quintiere 


Abstract 


An  analytical  model  was  developed  to  describe  the  processes  involved  in  the  burning  and  extinction 
of  solid  materials.  Included  are  flame  heat  transfer,  charring,  transient  conduction,  and  water 
application.  The  model  qualitatively  describes  the  burning  rate  of  both  charring  and  thermoplastic- 
like solids.  It  illustrates  how  the  steady-state  heat  of  gasification  can  be  derived  from  peak  burning 
rate  test  data  taken  as  a function  of  irradiance.  Experimental  data  are  shown  to  support  this 
derivation.  The  model,  in  conjunction  with  a critical  flame  temperature,  is  used  to  describe 
suppression  and  extinction  by  water. 

Keywords:  Burning  rate,  extinction,  heat  of  gasification,  fire  models. 


1.  Introduction 

The  degradation  of  a solid  in  burning  is  a complex  process.  It  involves  flame  heat  transfer, 
evaporation  or  pyrolysis,  possibly  charring,  and  transient  thermal  effects.  Although  specific  models 
have  been  developed,  they  are  limited  to  classes  of  materials  and  do  not  apply  in  general.  Test 
apparatuses,  such  as  the  Cone  Calorimeter,  provide  a means  for  dynamically  measuring  the  mass  loss 
and  energy  release  of  solid  materials,  but  the  interpretation  of  the  data  is  limited  by  the  lack  of  a 
simple  burning  rate  model.  The  purpose  of  this  presentation  is  to  describe  a simple  model  that 
represents  the  significant  effects  in  the  burning  of  a solid,  and  to  use  that  as  a potential  means  for 
analyzing  burning  rate  data.  The  intent  is  to  lay  out  the  governing  equations  so  that  some  limited 
quantitative  insight  might  be  derived.  The  literature  has  not  been  reviewed  for  other  charring  models 
since  depth  of  specific  detail  is  not  the  objective  here.  Instead  a simple  but  complete  model  for  the 
principal  factors  of  solid  combustion  are  sought.  At  some  later  time,  solutions  to  the  governing 
equations  will  be  conducted,  and  their  relationship  to  data  and  other  models  in  the  literature  will  be 
examined.  Ultimately  we  would  like  to  have  a model  that  is  universally  suited  for  describing  the  bum 
rate  of  any  complex  material  in  term  of  experimentally  derived  parameters. 


2.  Model 
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A one  dimensional  steady  "stagnant-film"  approach  is  used  to  describe  the  gas-phase  combustion 
processes.  The  "stagnant-film"  implies  that  the  fluid  boundary  layer  will  be  at  rest,  but  its  thickness 
will  be  determined  from  known  heat  transfer  relationships.  Terms  for  flame  radiation  and  external 
radiation  are  included.  It  is  assumed  that  variables  respond  quickly  in  the  gas  phase  so  that  it  is 
quasi-steady  with  respect  to  the  solid  phase.  In  the  solid  phase,  charring,  vaporization  and  transient 
conduction  effects  are  considered.  As  vaporization  of  fuel  occurs  in  the  solid  it  instantaneously  is 
transported  to  the  surface  with  no  resistance  through  the  char  layer.  A global  analj^is  will  be  used 
to  describe  the  solid  phase  processes  which  will  be  linked  to  the  gas-phase  by  the  surface  boundary 
condition.  The  stagnant  film  model  is  available  in  standard  combustion  texts  (A.M.  Kanury, 
Introduction  to  Combustion  Phenomena.  Gordon  and  Breach,  1977,  pp  167-178),  but  is  repeated  here 
for  completeness.  First  the  solid  phase  analysis  will  be  considered. 

Figure  1 illustrates  the  model  for  decomposition  in  the  solid.  Vaporization  is  assumed  to  occur  at 
a plane  behind  the  char  zone.  Figure  2 decomposes  this  model  into  control  volume  regions,  and 
displays  the  heat  and  mass  transfer  processes  for  each  region.  Measurements  are  based  on  the  mass 
loss  rate  of  the  solid  during  burning,  so  we  seek  to  represent  the  processes  in  terms  of  this  quantity. 
Moreover,  the  heat  of  combustion  or  specie  yields  derived  from  measurements  are  given  with  respect 
to  the  measured  mass  loss  rate.  Thus,  one  may  interpret  these  properties  as  effective  values  for  the 
fuel  gaseous  mixture  that  leaves  the  solid.  For  example,  if  water  vapor  is  driven  off  in  addition  to 
fuel  pyrolyzate,  the  heat  of  combustion  would  represent  that  of  this  mixture. 

In  the  following  sections  a detailed  derivation  of  the  solid  and  gas  phase  equations  will  be  presented. 
For  the  most  part,  the  control  volume  conservation  laws  for  mass  and  energy  will  utilized.  As  a check 
and  to  offer  another  perspective,  a differential  formulation  of  the  decomposing  solid  will  be 
examined.  The  results  will  be  shown  consistent  with  the  control  volume  formulation. 

2.1  Solid  Phase  Model 

The  conservation  of  mass  and  energy  will  be  applied  to  the  regions  depicted  in  Figures  2a-2c.  These 
consist  of  the  (a)  the  char,  (b)  the  vaporization  interface,  and  (c)  the  virgin  solid.  For  mass 
conservation  with  respect  to  a control  volume,  the  rate  of  change  of  mass  within  the  control  volume 
plus  the  net  rate  of  mass  flow  out  of  the  contol  volume  is  equal  to  zero.  The  conservation  of  energy 
for  a control  volume  states  that  the  rate  of  change  of  internal  energy  within  the  contol  volume  plus 
the  net  rate  of  enthalpy  out  is  equal  to  the  net  rate  of  heat  addition.  Work  is  neglected,  except  for 
flow  work,  and  combustion  energy  can  be  treated  as  an  effective  heat  addition.  Constant  properties 
will  be  assumed  for  each  distinct  media,  i.e.  char  and  virgin  material;  and  the  internal  energy  is 
referenced  at  the  initial  temperature. 

Conservation  of  mass  for  the  entire  solid 


The  conservation  of  mass  for  the  solid  in  Figure  1 is  given  by 


1 dm  • t!  r\ 

+ iw"  = 0 G) 

A dt 
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where 


— is  the  mass  rate  of  change  of  the  solid, 


dt 


A is  the  surface  area, 

and  ih"  is  the  mass  flow  rate  of  the  gasified  products  per  unit  area. 


Hence,  rii"  is  the  mass  loss  rate  per  unit  surface  area,  or  which  is  more  commomly  referred  to  as  the 
burning  rate  (per  unit  area).  Eq.  (1)  is  not  essential  to  the  governing  equations  which  follow,  but  is 
only  presented  to  show  the  relationship  between  the  rate  of  fuel  gasses  which  flow  out  of  the  system 
and  the  rate  of  mass  loss  of  the  solid. 

Conservation  of  mass  for  the  char 

By  considering  the  conservation  of  mass  for  the  char  layer  in  Figure  2a,  and  assuming  that  the 
vaporized  fuel  instantaneously  leaves  the  solid,  it  follows  that 


he 


(2) 


where  is  the  char  density,  and  rtij."  is  the  char  mass  flow  rate  per  unit  area  from  the  vaporization 
plane.  The  accumulation  of  the  fuel  gases  in  the  char  matrix  is  taken  to  be  negligible,  so  that  the 
rate  of  gaseous  fuel  entering  from  the  vaporization  interface  is  equal  to  the  rate  of  fuel  gases  which 
leave  the  solid  to  be  be  burned. 

Conservation  of  mass  for  the  virgin  solid 

Consequently,  a conservation  of  mass  on  the  virgin  fuel  element  in  Figure  2c  yields 


(3a) 


pdy^  + pv  = 0 


where  p is  the  density  of  the  virgin  solid,  and  where  the  control  volume  surface  moves  at  the  speed 

d6  db 

of  the  vaporization  plane,  v = — - = - — - . It  is  assumed  that  the  thickness  of  the  solid  stays 

dt  dt 

constant  as  decomposition  occurs. 

Conservation  of  mass  for  the  vaporization  interface 
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A mass  balance  at  the  vaporization  plane  together  with  Eq.(2),  and  assuming  is  constant,  gives 


pv  = m”  + m"  = pw" 
where  p = p/(p-p> 


(3b) 


This  equation  states  that  the  rate  of  mass  from  the  virgin  solid  into  the  vaporization  interface  is  equal 
to  the  mass  flow  rate  of  char  and  fuel  gases  leaving  the  vaporization  interface. 

The  conservation  of  energy  for  each  region  in  Figure  2 is  now  considered. 

Conservation  of  energy  for  the  char 

For  the  char  region  (2a), 


(4) 


- aT’/  - q 


. // 


where 


Uc 

is  the  internal  energy  of  the  char  layer. 

h 

is  the  specific  heat  of  the  vaporized  fuel  gas. 

Cc 

is  the  char  specific  heat. 

-Ik 

is  the  flame  convective  heat  flux. 

r 

ctT 

is  the  flame  radiative  heat  flux, 

is  the  external  radiative  flux  to  the  surface, 

is  the  surface  reradiative  heat  flux  (assumed  to  be  a blackbody). 

Aft 

^ V 

is  the  heat  flux  to  the  plane  of  vaporization. 

T 

is  the  vaporization  temperature. 

T 

is  the  surface  temperature. 

T 

is  the  reference  temperature. 

Conservation  of  energy  for  the  vaporization  interface 

The  conservation  of  energy  applied  to  the  vaporization  plane  of  Figure  2b  yields: 
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(5) 


/ \ A TJ  'll  'll 

(pv)A/r^  = 

where  AH^  is  the  heat  of  vaporization  (pyrolysis)  for  the  solid  at  temperature  T^, 
is  the  heat  loss  per  unit  area  to  the  virgin  solid. 

The  left  side  of  Eq.  (5)  relates  to  the  energy  required  to  change  the  virgin  solid  to  vapor  and  char, 
and  can  be  taken  as  a definition  of  AH^ 

Conservation  of  energy  for  the  virgin  solid 

An  energy  balance  on  the  virgin  solid  gives 


1 ^ 

A dt 


+ pvc(r,-r^  = q'^  - qH 


(6) 


where  is  the  heat  loss  per  unit  area  &om  the  back  of  the  solid, 

and  Uy  is  the  internal  energy  of  the  virgin  solid. 

Since  the  density  (p)  and  specific  heat  (c)  of  the  virgin  solid  can  be  considered  constants, 

U,  = pc4  f (.T-T,)dy,  CO 

«. 


Equations  (2)  - (7)  constitute  the  governing  equations  for  the  solid  phase.  It  is  not  entirely  obvious 
what  are  the  unknown  variables,  and  what  is  the  strategy  of  solving  for  them.  For  now  it  can  be 
noted  that  by  assuming  appropriate  proGles  for  the  temperature,  the  heat  fluxes  can  be  expressed  in 
terms  of  temperature  by  Fourier’s  Law.  Also  the  variables  involving  the  gas  phase  heat  transfer  need 
to  be  developed  firom  the  gas  phase  analysis  to  follow  in  the  next  section.  Now  a digression  is 
introduced  to  illustrate  another  approach  for  presenting  the  solid  phase  equations,  and  to  offer  a 
check  on  the  above  analysis. 

Differential  formulation 


Consider  pure  conductive  heat  transfer  into  the  virgin  solid  with  the  space  coordinate  (x),  fixed  to 
the  moving  vaporization  plane.  This  coordinate  system  is  shown  in  Figure  1,  and  is  introduced  to 
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avoid  confusion  with  the  coordinate  system  used  above.  If  is  the  initial  fixed  reference  system, 


X = 


- fv(x)dT 


0 


(8) 


where  v is  the  velocity  of  the  vaporization  plane.  The  conduction  equation  in  the  fixed  frame  of 
reference  is 


which  transforms  as 


(9) 


(10) 


from  Eq.  (8).  Hence  in  the  moving  frame  of  reference, 


(11) 


with  the  conditions: 
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and 


f=0,  r=r^  the  initial  temperature. 


Case  1.  Non-charrinp  steady  burning 

Let  us  consider  the  ideal  case  of  a non-charring  material  undergoing  steady  burning.  If  steady 
conditions  prevail  in  the  moving  system,  i.e.,  the  temperature  field  is  not  changing  in  the  virgin  solid 
relative  to  the  moving  vaporization  plane,  and  the  back  face  conditions  are  negligible,  i.e.,  a very  thick 
solid,  then  Eq.  (11)  becomes 


From  Eq.  (3)  and  since  v = - — - 


dt 


(14) 


with  conditions  from  Eq.  (5)  and  (12) 


(15a) 


X = 0,  r = 


(15b) 


X 


(15c) 


Using  Eqns  (15b)  and  (15c) 


= e 


(16) 


T-T^ 

Ty-T, 


and  from  Eq.  (15a) 


^li- 
nt - 


.// 

Qy 


AH^  * ciT^-T;^ 


(17) 


The  denominator  of  Eq.  (17)  is  commonly  referred  to  as  the  steady  state  heat  of  gasification  (L.). 


= AH^  * c(T^-TJ 


(18) 


Case  2.  Transient  charring 

In  general,  other  terms  which  will  be  considered  below  will  affect  the  mass  loss  rate,  rii".  If  the 
process  is  not  steady,  we  can  consider  Eq.  (11)  by  integrating  each  item  over  6^: 


= (-4i)  - i-ii!) 


-cm"  (T,  - r.) 


© V / 0 

K 

= pc—  r T dx  + cm" 

dt  J 
0 

0 
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Also  for  the  charring  case  since  at  the  vaporization  plane 
db  db 

V ^ = — - , and  from  Eq.  (3b)  it  follows  that  ih"  should  be  replaced  by  jSih"  in  the  above. 

dt  dt 

Substituting  in  the  integrated  form  of  Eq.  (11)  gives 


pc^j  (T-T,)dx  * ^H,"c(T,-T,)  - 

o 

. //  . // 

= 9*  - 9* 

which  is  identical  to  Eq.  (6).  This  demonstrates  that  the  conservation  of  energy  for  the  virgin  solid 
is  consistent  with  this  differential  formulation. 

Solution  strategy 

In  review,  six  independent  equations  have  been  presented  consisting  of  conservation  of  mass  and 
energy  for  the  char,  vaporizing  interface,  and  virgin  solid.  These  can  be  combined  so  that  only  two 
energy  equations  (4  and  6)  are  considered  with  the  unknown  variables:  char  depth,  surface 
temperature,  and  a variable  to  be  introduced  which  represents  the  thermal  effects  in  the  virgin  solid. 
The  thickness  of  the  virgin  solid  and  the  mass  burning  rate  can  be  related  to  the  char  depth,  and  the 
heat  fluxes  can  be  related  to  temperature  profiles  to  be  assumed  in  terms  of  the  thicknesses  of  the 
char  and  virgin  solid.  A third  equation  will  be  found  from  the  gas  phase  flame  analysis  which  will  give 
a relationship  for  the  burning  rate.  The  flame  convective  and  radiative  heat  fluxes  will  be  expressed 
in  terms  of  the  current  variables  from  the  gas  phase  flame  analysis.  Below,  this  section  will  be 
concluded  with  a derivation  of  the  convective  flame  heat  transfer  based  on  the  above  conservation 
equations  which  will  serve  as  a boundary  condition  to  the  gas  phase  analysis  developed  in  the  next 
section.  This  couples  the  solid  and  gas  phase  analyses. 

From  Eqns.  (5)  and  (6) 


A 

pci  / (r-r.)&  * pi>>"c(r,-r.)  = q'J 

or  ' (19) 

= i"  - 4i'  - pc|  / (r-rj  dx. 
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This  is  a departure  from  the  steady  state  results  given  by  Eq.  (17).  It  applies  to  the  non-charring  case 
if  is  regarded  as  the  net  surface  heat  flux,  and  /S  = 1. 

The  more  general  charring  case  is  considered  by  eliminating  from  Eq.  (19)  by  using  Eq.  (4). 
Furthermore,  the  internal  energy  of  the  char  is  represented  as  follows: 


A 


= / Pc^C^S 


where  y^  is  measured  from  the  charring  surface.  The  internal  energy  per  unit 


mass,  Up,  is  represented  as  Cp(T-Tq)  where  Cp  is  the  char  specific  heat.  By  Eq.  (2)  it  follows  that 


idU^ 
A 


dU.  A t 


(20) 


Combining  all  of  the  energy  equations  for  the  solid,  or  from  Eqns.  (4),  (19)  and  (20) 


. •> 

/ (r-r.)  dy,  - 

e 

V c 

-4;  - / (T-T.)dy. 

- pM\  (r,-rj 


(21a) 
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Alternatively, 


Pm%  = 


dT,^ 

[ ‘ dy) 


+ 

+ 

a" 

Hext,r 

ar/ 

fymc 

flame 

external 

reradiation 

radiation 

radiation 

(noncharring) 

-r.‘)  - 

d ) 

[T-T,)dy. 

surface  reradiation  energy  storage 
due  to  charring  due  to  charring 


energy  flow  through  char  back  face  heat  loss 


(21b) 


-pc 


virgin  solid  energy  storage 


(The  above  labels  are  qualitative  descriptions  of  the  terms.) 
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Equation  (21b)  gives  the  thermal  boundary  condition  for  the  gas  phase  analysis,  i.e., 


dT 

- k(—^)  = Pm"  Sf 
* dy 


pc  j-  ) (X-T,)dy,  * qH  * j(T  - T^dy, 

where  a s i,  + 

‘ Pm" 

^ m\(r,  - rj  - p,M"c,  (7;-r,)/p  ^ ot^  (22) 

pm" 

. //  . // 

_ <lf^^  Qext^ 

Pm" 


The  form  in  Eq.  (22)  constitutes  a boundary  condition  for  the  gas  phase  problem  to  follow. 


2.2  Flame  Model 


The  gas-phase  one  dimensional  diffusion  flame  is  solved  by  conventional  means  (e.g.,  Kanury  [7]). 
The  conservation  of  mass,  energy  and  species  are  considered  for  a vaporizing  solid.  A stagnant  film 
of  thickness  6g  represents  the  boundary  layer  thickness  or  domain  of  the  combustion  region.  Only 
one  dimensional  diffusive  effects  are  considered  with  all  of  the  convective  effects  due  to  transverse 
fluid  motion  accounted  for  by  deriving  6g  from  known  convective  heat  transfer  results.  The  steady 
equations  are  listed  below  and  relate  to  the  diagram  in  Figure  3: 


12 


Conservation  of  mass 


m‘'  = constant  for  all  y , 


(23) 


Conservation  of  energy 


nt  c. 


fh 

dy 


= k. 


d^T 


dy'^ 


i _ 


/// 


m/'AHfl-x;) 


(24) 


Conservation  of  oxygen 


„ dy„ 

m"  — = p.  C 2£ 

dy  ‘ dy^ 


+ m. 


/// 


(25) 


All  of  the  fluid  properties  are  assumed  constant.  D is  the  diffusion  coefOcient  and  the  Lewis  number, 
,gCgD  is  assumed  equal  to  one.  Stoichiometry  gives 


. ///  . /// 
= rm^ 


(26) 


where  r is  the  stoichiometric  oxygen  to  fuel  mass  ratio,  in"  and  rfa"jj  are  the  respective  rates  of 
generation  per  unit  volume  for  fuel  and  oxygen. 

Boundary  conditions 
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The  boundary  conditions  are  given  as 


y-\,T-T. 

= JW- 

y = 0,  7-  = r. 


dY. 


ox 


dy 


= 0, 


dT 

and  = m"pa!  Eq.  (22). 

^ dy 


(27) 


Burning  rate  solution 
By  introducing  a new  variable, 


b = T * 


A^,(l  - 


(28) 


based  on  adding  Eqns.  (24)  and  (25)  in  order  to  eliminate  the  chemical  source  terms  it  can  be  shown 
that 


w"  = — In  (1  +B) 


(29) 


where 


_ ^ (r.  - r,) 

P9! 


(30) 


and 


has  been  replaced  by  h,  a convective  heat  transfer  coefficient. 
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Flame  heat  flux 


By  introducing 


$ 


in  c. 


(31) 


it  can  be  shown  from  Eq.  (22)  that  at  y = 0 


dT, 

K— 

^ dy 


? 1 

h 

A^,(l-X.) 


ox,* 


+ cXT^  - T,) 


(32) 


which  is  the  convective  heat  flux  to  the  surface.  The  quantity  5/(^^-l)  is  called  the  "blocking 
factor"  which  effectively  reduces  h from  its  pure  heat  transfer  value  as  the  mass  transfer,  ih", 
increases. 

Closure  on  solution  strategy 


Eqns.  (29)  and  (32)  provide  the  two  needed  additional  relationships  that  allow  a solution  to  the 
complete  set  of  equations.  A solution  of  the  transient  equations  will  not  be  attempted  in  this 
presentation,  but  instead  some  qualitative  observations  will  be  made,  and  an  extinction  criterion  will 
be  examined.  In  a future  investigation  it  is  intended  that  these  equations  will  be  solved,  and  the 
extinction  criterion  will  be  used  to  obtain  a full  simulation  of  the  burning  rate  for  a typical  solid. 


3.  Qualitative  Results 

Although  it  is  possible  to  pursue  approximate  solutions  by  assuming  temperature  profiles  for  the  solid, 
this  will  not  be  done.  Instead,  qualitative  analyses  will  be  explored  to  describe  the  general  burning 
behavior  of  solids.  Initially,  as  the  solid  reaches  its  vaporization  temperature  Ty  it  behaves  like  a 
thermoplastic  or  non-charring  material  before  a significant  char  layer  is  developed.  To  simplify  the 
interpretation  it  is  useful  to  initially  consider  the  blocking  factor  and  the  Oame  radiant  heat  flux  to 
be  constants.  The  former  will  decrease  with  ih"  while  (Jf , is  likely  to  increase  with  ih"  due  to  an 
increase  in  the  flame  thickness.  Both  of  these  effects  can  oe  used  to  temper  the  descriptions  below. 

Figure  4 depicts  the  behavior  of  the  burning  rate  over  time  as  vaporization  ensues  after  time  zero. 
For  a non-charring  material,  the  surface  approximately  remains  at  Ty  and  only  the  energy  storage 
term  in  the  material  affects  the  bum  rate  provided  the  solid  is  thick  so  no  significant  back-face  heat 
loss  occurs.  As  time  increases,  the  energy  storage  term  approaches  zero  (see  Eqs  (17)  and  (19))  as 
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a steady-state  is  sought.  However,  once  thermal  effects  are  felt  at  the  back-face,  the  bum  rate  is 
disturbed.  For  a back-face  substrate  with  a lower  conductivity  than  the  burning  solid,  an  increase  in 
the  bum  rate  will  occur  until  bum-out  of  the  solid.  The  opposite  effect  will  result  for  a higher 
conductivity  substrate. 

If  charring  occurs  the  behavior  is  more  complex,  but  describable  by  Eq.  (21b).  Initially,  the  burning 
begins  like  the  non-charring  case.  But  as  the  char  layer  increases,  increases  and  the  charring  terms 
in  Eq.  (21b)  all  increase  to  reduce  the  burning  rate.  This  is  shown  in  Figure  4,  again  with  a 
corresponding  back-face  heat  loss  effect.  At  the  maximum  for  the  charring  curve  and  for  the  peak 
of  the  non-charring  curve  before  back-face  effects  occur,  it  can  be  argued  that  the  following  terms 
are  small: 


8 4-8 

d t 

pc  — / (T-T^)  dy^  « 0 (since  we  assume  the  peak  is  near  steady-state) 

dt  { 

8, 


o (r/  -T*)  * 0 (since  before  charring  effects  may  increase  the  surface 

significantly  above  its  ignition  temperature) 


P c 

c 


(T  -T^  dy^  - 0 (since  the  char  layer  is  thin) 


cfX^  - 7’^)/p  = 0 (if  non-charring  or  if  p^  is  small) 


m"c^(T^-T^)  « 0 (since  T,  « T^) 


q'l  e 0 (since  the  solid  is  assumed  thermally  thick  which  would 
hold  in  the  early  stage  of  burning  for  materials) 


Admittedly,  the  above  arguments  are  not  subtantiated,  and  must  require  analysis  of  at  least  the 
solutions  for  the  burning  rate  problem.  However  they  are  accepted  as  approximations  for  now. 
Therefore,  these  peak  burning  rate  conditions,  Eq.  (21b)  can  be  approximated  as 
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m. 


e ] 

* f 

ye  l ) 

c 1 

g 

[n»,.  iffcd  - xp/r-c,(r.  -T^]*  4;!,^  q'L,  - ° 


(33) 


where  jSLg  can  be  regarded  as  the  steady  state  heat  of  gasification  for  a charring  materials  in  general. 

If  in  a test  apparatus  the  convective  and  radiative  flame  heat  fluxes  do  not  vary  greatly  over  a range 
of  peak  mass  loss  rates,  then  the  variation  of  with  ^ ^ linear.  The  slope  is  IZ/SL.  and 

provides  an  approximate  method  for  determining  /SLg.  Moreover,  if  the  flame  heat  fluxes  could  be 
determined,  (33)  would  provide  a means  for  estimating  the  peak  bum  rate  per  unit  area  for,  in 
principle,  all  fire  conditions.  Thus,  given  jSLg  for  the  material,  and  a knowledge  of  the  fire  heat 
fluxes,  the  burning  rate  can  be  estimated.  This  procedure  provides  a framework  for  utilizing  small 
scale  data  such  as  that  developed  from  the  Cone  Calorimeter  and  similar  devices. 


3.1  Examination  of  Experimental  Results 

It  is  generally  found  that  small  scale  test  data  for  the  mass  or  energy  release  rate  of  materials  is  an 
approximate  linear  function  of  irradiance  ^ Tewarson  [1]  has  shown  this  to  be  the  case,  and  has 
utilized  heating  in  nitrogen  to  derive  Lg  from  Eq.  (33),  In  the  case  of  nitrogen,  the  flame  heat 
transfer  terms  are  non-existent.  Jackson  [2]  performed  similar  experiments  in  a nitrogen  atmosphere 
and  investigated  effects  over  time.  He  defined  a bulk  heat  of  gasification  of  the  solid  as  the  total  net 
heat  added  divided  by  the  instantaneous  mass  loss  rate.  In  terms  of  Eq.  (21),  this  would  correspond 
in  a nitrogen  (non-flaming)  system  as 


A 

(34) 

- p.Pc.(n-7’«)/p 


Based  on  our  previous  discussion,  the  bulk  heat  of  gasification  for  a non-charring  material  would 
asymptotically  approach  Lg,  the  steady  value  as  the  energy  storage  term  goes  to  zero.  Jackson, 
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indeed,  finds  this  to  occur  for  PMMA  measurements  fi-om  which  yield  (Lg)buik  as  a function  of  time, 
but  the  asymptote  depended  on  the  irradiance  level.  This  dependence  could  have  been  a result  of 
inaccuracies  in  his  estimation  of  the  back-face  heat  loss  terms  as  suggested  by  the  fact  that  the  effect 
on  Lg  was  diminished  as  the  irradiance  increased.  Figure  5 shows  Jackson’s  results  for  PMMA  in 
terms  of  peak  mass  loss  rate  plotted  for  various  irradiation  levels.  The  slope  of  those  data  yield  L. 
= 1.89  kJ/g  as  suggested  by  using  Eq.  (33).  This  is  compared  to  a time  asymptotic  value  at  40  kW/m^ 
irradiance  of  1.65  kJ/g.  TTiis  compares  to  a steady-state  value  by  Tewarson  of  1.63  kJ/g  [1].  Figure 
6 shows  peak  values  of  m"  plotted  for  oak  as  taken  by  Jackson  in  nitrogen.  From  the  slope  of  these 
data,  Lg  (or  more  precisely  /3Lg)  for  oak  is  determined  as  4.0  kJ/g. 

An  example  of  data  taken  under  flaming  conditions  is  presented  in  Figure  7 for  particle  board.  These 
samples  were  tested  vertically  at  a nominal  size  of  0.3  m high.  The  data  come  from  Walton  and 
Twilley  [3]  for  irradiances  of  25,  50  and  74  kW/m^,  and  for  0 irradiance  as  determined  by  Kulkami 
[4].  Also  shown  are  data  from  Tsantiadis  and  Ostman  [5]  for  a particle  board  material  that  may  not 
necessarily  be  identical  to  that  of  the  other  studies.  Whether  the  absolute  peak  or  an  average  peak 
value  for  ih"  is  used,  the  results  for  Lg  range  fi-om  4.4  to  4.7  kJ/g.  If  one  applies  this  technique  to 
a wide  range  of  data  for  materials  used  in  a series  of  room  fire  tests  performed  in  Sweden  [5],  the 
Lg  values  can  be  determined  and  are  given  in  Table  1.  It  is  not  possible  to  develop  an  explanation 
for  the  trends  of  these  data,  and  to  explain  them  in  terms  of  their  material  composition.  The  result 
for  Lg  can  simply  be  regarded  as  empirical  which  allow  the  prediction  of  the  change  in  mass  loss  rate 
with  me  change  in  surface  heat  flux.  However,  Eq.  (33)  gives  a physical  interpretation  to  Lg  making 
it  the  steady-state  heat  of  gasification  appropriate  to  burning  a liquid-like  fuel.  Although  this 
mechanism  of  pyrolysis  will  be  very  different  from  the  evaporation  of  a liquid,  the  model  being  used 
appears  to  give  a first-order  representation  of  peak  mass  loss  with  heat  flux  in  a range  of  conditions 
appropriate  to  common  building  fires. 


4.  Extinction 

To  close  this  simplified  view  of  burning,  extinction  is  considered  in  terms  of  a critical  flame 
temperature  criterion.  It  is  empirically  established  that  for  both  pre-mixed  flames  at  the  lower 
flammable  limit  and  for  diffusion  flames  at  extinction,  a minimum  flame  temperature  corresponds  to 
these  phenomena.  Typically  this  critical  temperature  is  roughly  1600  K for  many  C-H-O  fuels  burning 
in  air  (e.q.  >\^lliams  reports  measured  values  at  extinction  of  1500  ± 50K  [6]).  Based  on  this 
concept,  it  can  be  investigated  how  the  flame  temperature  in  the  stagnant  film  flame  model  can  be 
reduced  through  heat  losses.  Equation  (21)  can  be  used  as  a boundary  condition  to  establish  the  net 
heat  losses  from  the  flame  at  the  fuel  surface.  One  approach  is  to  return  to  the  stagnant  film  flame 
model  and  solve  for  the  temperature  at  the  flame  sheet. 

Flame  temperature 

Because  we  have  included  more  heat  transfer  processes  than  would  be  consistent  with  the  pure 
convective  stagnant  film  model  we  must  consider  the  relevance  of  these  fluxes  on  flame  temperature. 
A control  volume  approach  is  taken  as  shown  in  Figure  3.  The  heat  and  mass  transfer  processes  that 
directly  affect  the  flame  temperature  are  illustrated  in  Figure  8.  The  radiative  loss  from  the  flame 
is  accounted  for  by  an  effective  decrease  to  the  heat  of  combustion  --  yielding  less  chemical  energy 
available  to  produce  the  flame  temperature.  Therefore,  the  only  heat  loss  from  the  flame  (control 
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volume)  involves  the  convective  heat  flux  to  the  fuel  surface  --  this  constitutes  shown  in  Figure 
8.  The  other  heat  fluxes  to  the  surface,  and  aT^,  from  the  surface  are  considered  not  to 
interact  with  the  flame  system. 

By  conservation  of  mass  under  steady  state  conditions,  the  mass  flux  of  combustion  products  flowing 
out  of  the  system  is  equal  to  the  mass  flux  of  fuel  gases  and  the  mass  flux  of  air  flowing  in  to  the 
system: 


= m 


// 


. // 

+ m 


a * 


(35) 


and  by  conservation  of  energy  and  Eq.  (22), 

mUc^CTf-TJ  - - TJ  = AH,(1-X,)  - ih"  pSf'  (36) 


. // 

where  + — — by  Eq.  (22)  because  the  flame  heat  flux,  (^"fj  has  been  effectively  eliminated 

pm '' 

in  this  analysis  since  all  the  flame  radiative  loss  is  accounted  for  by  Xp 
Also  by  stoichiometry: 


. //  . //  r 

= m 


and 


ox," 


m"  (1 


(37) 


Then  solving  for  Tf  yields 


[AH,(1  -x^  - pa'iy^./r  ^ c,  (r.  - r.) 


(38) 


which  can  be  found  in  Kanury  [7]  for  the  case  of  pure  convection:  X'  = Lg,  T^  = T^  and  = 0. 
The  quantities  that  reduce  Tf  and  drive  the  flame  towards  extinction  are  readily  seen: 
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1.  increasing 

2.  decreasing 

3.  decreasing 

4.  increasing  energy  storage  due  to  charring  and  transient  effects  in  the  solid 

5.  increasing  Lg 

6.  decreasing  external  radiative  heat  flux 

Extinction  by  the  evaporation  of  water  can  also  be  accounted  for  by  adding  ih^  IVni"  to  the  right 
side  of  Eq.  (22)  if  ih^  is  the  mass  evaporation  rate  of  water  per  area  of  the  solid  surface  and  is 
the  heat  of  gasification  of  the  water  (2,6  kJ/g). 

As  an  example  of  how  to  apply  the  temperature  criterion  consider  some  simple  cases. 

Example  1.  Steady  burning  of  PMMA 

Take  the  fuel  to  be  PMMA  and  ignore  all  transient  and  radiative  effects.  We  might  expect  flame 
radiative  effects  to  be  small  near  extinction,  but  not  transient  and  reradiation  effects.  Hence  our 
result  will  be  illustrative,  but  approximate.  The  following  properties  are  used: 

= 1 J/gK 

= 20“C 

= 377 °C 

= 25.2  kJ/g 

= 1.63  kJ/g 

= 2 g oxygen/g  PMMA 

Using  Tf  = 1600  K,  we  can  estimate  the  Yqj,oo  to  cause  extinction  of  this  PMMA  flame.  Substituting 

into  Eq.  (38)  yields: 


(25.2  - 1.63)r  ./2  + (10-^)  (293  - 650) 


(10-^)  (1600  - 650)  = 


ox: 


1 " Yox,J^ 


Yox,.= 


1.376 

4.79 


= 0.117, 


or  extinction  occurs  at  a mole  fraction  of  approximately  10.5%  oxygen.  The  corresponding  burning 
rate  at  extinction  in  this  oxygen  atmosphere  can  be  found  from  Eq.  (33)  as  follows: 

Let  h = 14  W/m^K  (a  typical  value  corresponding  to  natural  convection)  and  assume  the  blocking 
factor  = 1. 


m 


// 

extinct 


+ 


- r^]  - a 


d/i, 
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where  now  we  include  the  re-radiation  loss. 


= |(y)  [0.1 17  (25.2)/ 2 + 10”^  (293  -650))  - 5.67x10-^1(650)^1/  1.63 


or 


m 


// 

extinct 


3.4  glm^s . 


The  heat  transfer  coefficient  used  above  was  selected  as  a plausible  value  that  would  lead  to  a typical 
value  for  ih"  found  at  extinction.  Thus,  under  the  condition  given  above,  the  PMMA  will  extinguish 
in  an  oxygen  environment  of  10.5%  with  a maximum  bum  rate  of  3.4  g/m^s  just  before  extinction. 

By  considering  all  the  terms  in  Eq.  (38)  we  can  determine  a relationship  among  the  "independent" 
parameters  that  must  hold  at  extinction.  For  example,  for  fixed  values  of  T^,  Y^jp^,  T^,  there  is  a 
critical  value  of  5E*  at  extinction  (Tf  = 1600K). 

Example  2.  Ouasi-steadv  suppression  of  PMMA 

As  another  sample  let  us  examine  the  steady  burning  of  PMMA  simultaneously  exposed  to  an 
irradiance  ((J  and  a steady  water  spray  (ih  ^).  From  Eqns.  (22)  and  (36), 


S£'  = i,  - 


7A  , . /fw  . // 

V ~ 9ext,r 


m 


// 


(39) 


Let  Yqjpoo  = 0.233,  T^  = 293  K and  the  previous  PMMA  properties.  Substituting  into  Eq.  (38) 
yields 


10-^  (1600  - 650) 


_ {[25.2 -(1.63  + 5.67x10-^^(650)^  + m^^(2.6)-4^^]} 


m 


extinct 


.//  _ 2.6m^^  10.1  - 

or 


(40) 
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From  the  steady  PMMA  (non-charring)  mass  loss  Equation  (21)  with  the  water  term  included  as 
on  the  right  side, 


'h} 


. // 


.// 

^extjr 


-ot;  - 


• H T 


where  the  blocking  factor  has  been  taken  as  one.  Here  we  will  take  h = 8 W/m^K  to  match  the  data 
for  no  water  application.  Ignoring  flame  radiation  loss,  yields 

m"(L63)  = (|)  [0.233  (25.2)/ 2+10-^ (293  -650)]  + 4^,  - 10.1  - 2.6in^' 

(41) 

otrh"  = 6.4  0.61  (-f-) 

m^-s 


with  in  kWIm^.  Eq.  (41)  gives  the  burning  rate  under  all  conditions.  Equating  this  to  Eq.  (40) 
gives  the  relationship  holding  at  extinction: 


0.20in^'  + 0.77  - 0.076  = 6.4  + 0.61  qZ,r  " 1*6 

or  q^,r  = 2.6  - 8.2  (— ) 


(42) 


Eq.  (42)  gives  the  critical  irradiance  for  a given  water  spray  rate  at  extinction.  The  results  of  Eq.  (41) 
and  (42)  are  in  qualitative  agreement  with  the  experimental  data  of  Magee  and  Reitz  [8]  as  shown 
in  Figure  9 with  the  critical  irradiance  at  extinction  underestimated  by  the  theory.  Nevertheless,  these 
results  show  the  versatility  and  level  of  accuracy  of  the  simple  burning  rate  model  presented. 


5.  Conclusions 

A simple  burning  rate  model  has  been  developed  for  a solid  which  is  based  on  a stagnant  film 
diffusion  flame  and  a control  volume  formulation  for  the  solid.  External  irradiance,  charring  and 
other  transient  effects,  back-face  heat  loss  and  water  application  are  phenomena  included.  The 
model  is  shown  to  yield  results  that  could  have  acceptable  accuracy  for  making  estimates  of  burning 
rate  for  materials.  It  is  shown  that  the  steady-state  heat  of  gasification  could  be  estimated  firom  small- 
scale  test  data  by  examining  the  slope  of  the  peak  burning  rate  versus  irradiance  for  both  charring 
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and  non-charring  materials.  The  terms  in  the  model  qualitatively  describe  the  burning  rate  seen  in 
test  data  for  both  charring  and  non-charring  materials.  An  extinction  model  based  on  a critical  flame 
temperature  yields  results  which  are  qualitatively  correct  with  accuracy  that  may  be  sufficient  to 
predict  first-order  effects. 
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Nomenclature 


A - area 

B -Eq.  (30) 

c - specific  test 

D - diffusion  coefficient 

h - convective  heat  transfer  coefficient 

k - thermal  conductivity 

Lg  - heat  of  gasification 

- heat  of  gasification  for  water 
^ - Eq.  (22) 

T - Eq.  (36) 
m - mass 
q - heat  transfer 

r - stoichiometric  oxygen  to  fuel  mass  ratio 
t - time 

T - temperature 
u - internal  energy  per  unit  mass 
U - internal  energy 

V - velocity 

X - coordinate 

- radiative  fraction  of  chemical  energy 
y - coordinate 

Yqj-  oxygen  mass  fraction 
jS  - density  ratio,  Eq.  (3b) 
h - thickness 
AHp-  heat  of  combustion 
AHy-  heat  of  vaporization 
p - density 
{ - Eq.  (31) 

Subscripts 

a - air 
b - backface 
c - char 
ext  - external 
f - flame 
g - gas 
k - conduction 
o - reference  state 
p - products 
s - solid,  surface 

V - vaporization,  virgin 
w - water 

00  - ambient 
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Superscripts 


(•)  - per  unit  time 
( )"  - per  unit  area 
( )*•  - per  unit  volume 
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Table  1 


Heat  of  GasiHcation  Values  at  Peak  Burning  Conditions 


Material 

(u^g) 

Data  Source 

PMMA 

1.89 

Jackson  [2] 

PMMA 

1.63 

Tewarson  [1] 

Oak 

4.0 

Jackson  [2] 

Particle  Board 

5.4 

[3],  [4],  [5] 

Insulating  Fiber  Board 

4.2 

Ostman  [5] 

Medium  Density  Fiber 

Board 

4.8 

tf  fl 

Spruce  panel 

6.3 

ff  fl 

Melamine-faced  Particle 

Board 

4.8 

tl  fl 

Paper  Wallcovering  on 

Gypsum  Board 

4.8 

fl  tl 

Plastic  Wallcovering  on 
Gypsum  Board 

3.7 

N n 

Paper  Wallcovering  on 

Particle  Board 

6.5 

II  n 

Textile  Wallcovering  on 
Gypsum  Board 

1.5 

fl  tl 

Textile  Wallcovering  on 
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Figure  1.  Solid  decomposition  model 
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Figure  2b.  Vaporization  process 
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Figure  2c.  Virgin  material  heat  transfer  processes 


29 


E 

o 


30 


Figure  3.  One  dimensional  diffusion  flame 
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Figure  4.  Qualitative  description  of  burning  rate  for  solids 
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Figure  5.  Peak  rates  of  mass  loss  versus  irradiance  for  PMMA  from  Jackson  [2] 
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Figure  6.  Peak  mass  loss  rate  for  oak  from  Jackson  [2] 
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Figure  7.  Peak  mass  loss  rate  for  particle  board  [3,4,5] 
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Figure  8,  Energy  and  mass  loss  streams  for  diffusion  flame 
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Figure  9.  Burning  rate  of  vertical  PMMA  slabs  versus  external  radiant 
flux  for  various  water  annlination  raw 
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